Dissipation effects in spin-Hall transport of electrons and holes 
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We investigate the spin-Hall effect of both electrons and holes in semiconductors using the Kubo 
formula in the correct zero-frequency limit taking into account the finite momentum relaxation 
time of carriers in real semiconductors. This approach allows to analyze the range of validity of 
recent theoretical findings. In particular, the spin-Hall conductivity vanishes for vanishing spin- 
orbit coupling if the correct zero- frequency limit is performed. 



I. INTRODUCTION 

In the recent years, an increasing interest in spin- 
dependent phenomena in semiconductors has developed, 
mostly in the field of spin electronics, which has by now 
become a major branch of solid state research 1-3 . On 
of the most investigated issues in this field is the influ- 
ence of spin-orbit coupling on various transport proper- 
ties of both electrons and holes. Many of these studies 
were inspired by the paradigmatic proposal of a spin field- 
effect transistor due to Datta and Das 4 ; for recent work 
in this direction see e.g. Refs. 5-13 . Most recently, in- 
teresting theoretical studies on the spin-Hall effect have 
been performed 14-16 . This effect amounts in a spin cur- 
rent (as opposed to a charge current) driven by an electric 
field perpendicular to it. In the present letter we reexam- 
ine these findings using the Kubo formula with full fre- 
quency dependence, treating both the case of electrons 15 
and holes 14,16 , and analyze the range of validity of pre- 
vious theoretical results obtained for the case of direct 
current. Here it is crucial to perform the correct zero- 
frequency limit taking into account an imaginary part of 
the frequencies occurring the the Kubo formula. 

The notion of the spin-Hall effect in systems of itin- 
erant spinful carriers was considered first by Dyakonov 
and Perel 17 in the early seventies, and, independently, in 
a more recent paper by Hirsch 18 . In these studies the pre- 
dicted spin-Hall effect is due to spin-orbit effects influenc- 
ing scattering processes upon static impurities. Following 
the terminology used in 15 ' 16 this is referred to as the ex- 
trinsic spin-Hall effect since it depends on impurity scat- 
tering. This is in contrast to the intrinsic spin-Hall effect 
predicted very recently in Refs. 1 which is entirely due 
to spin-orbit coupling terms in the single-particle carrier 
Hamiltonian and independent of any scattering process. 
As we shall see below, this distinction between intrinsic 
and extrinsic effects becomes ambiguous in the limit of 
weak spin-orbit coupling when life time effects of carrier 
quasiparticles have to be taken into account. 

Yet another type of spin-Hall effect was studied re- 
cently by Meier and Loss 19 in a two-dimensional Heisen- 
berg model consisting of isolated spins, in contrast to the 
itinerant-carrier systems mentioned before. 

Spin-orbit coupling also induces off-diagonal compo- 
nents of the conductivity tensor for charge transport. 



An important example is the anomalous Hall effect as 
it occurs in semiconductors in the presence of magnetic 
impurities 20 . Here, as in the case of the aforementioned 
spin-Hall effect, the off-diagonal elements of the charge 
conductivity tensor are the same in magnitude but differ 
in sign. Therefore, this antisymmetric conductivity ten- 
sor has the same components in all orthogonal coordinate 
systems, and in this sense the transport properties are 
isotropic. This is different from charge transport of elec- 
trons in quantum wells as investigated recently Ref. 12 . 
In such systems, the presence of spin-orbit coupling of 
both the Rashba 21 and the Dresselhaus 22 type leads to 
anisotropic dispersion relations and Fermi contours. This 
feature leads to symmetric off-diagonal elements in the 
conductivity tensor and therefore to preferred eigendirec- 
tions for charge transport 12,11 . This predicted effect of- 
fers a possibilty to detect spin-orbit coupling by measur- 
ing diffussive spin-unpolarized charge currents in a Hall- 
type geometry, which should be a comparatively simple 
experimental task. 

This paper is organized as follows. In section II we 
summarize elementary properties of linear-response the- 
ory as given by the Kubo formula. This technique is 
applied then in section III to spin-Hall transport of elec- 
trons in a quantum well in the presence of Rashba spin- 
orbit coupling. In section IV we investigate the case of 
bulk holes described by the Luttinger Hamiltonian in the 
spherical approximation. We end with conclusions in sec- 
tion V. 



II. KUBO FORMULA AND ZERO-FREQUENCY 
LIMIT 

Our present study of spin-Hall effect of electrons and 
holes in semiconductors is based on the usual Kubo for- 
mula with full frequency dependence for a spatially ho- 
mogeneous electric field 23 , 
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where we have assumed zero temperature T — and non- 
interacting carriers, which allows to formulate the two- 
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body Green's function entering the conductivity Kubo 
formula in terms of single-particle operators. A is the 
volume of the system, e is the elementary charge, and 
f(Sfj,(k)) is the T = Fermi distribution function for en- 
ergy £^,(k) at wave vector k in a dispersion branch labeled 
by [i. The velocity operators are given by v = i[H,r\/U 
where r is the position operator, and TL is the single- 
particle Hamiltonian not including the external electric 
field. The spin-current operator (in the Dirac picture) for 
spin moment polarized along the z-direction and flowing 
in the x-direction is given by 



with eigenstates 
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where s is the spin operator. The right hand side of 
Eq. (1) has to be understood in the limit of vanishing 
imaginary part rj > in the frequency argument. This 
imaginary part in the frequency reflects the fact that the 
external electric field is assumed to be switched on adi- 
abatically starting from the infinite past of the system, 
and it also ensures causality properties of the retarded 
Green's function occurring in Eq. (1). In general, and as 
we will discuss in more detail below, the limiting process 
rj — * does not commute with other limits, and, in par- 
ticular, the dc-limit u — > has to be taken with care 23 . 
In the presence of random impurity scattering, the re- 
tarded two-body Green's function in Eq. (1) will generi- 
cally have a frequency argument with positive imaginary 
part 23 . In this case the limit r\ — > is unproblematic, and 
the imaginary part of the frequency argument is just due 
to impurity scattering and/or other (many-body) effects. 
Generically, and as we will discuss in more detail below, 
the imaginary part r\ > corresponds to a finite carrier 
quasiparticle lifetime. 
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where xW = ay g(—k y + ik x ). By a straightforward cal- 
culation one obtains for the spin-Hall conductivity 
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are the Fermi momenta on the two dispersion branches 
for positive Fermi energy Ef > 0. In the presence of scat- 
tering on static random impurities, the imaginary part 
7] > in the frequency argument is given, to lowest order 
in the Rashba coefficient and the impurity potential, by 
the inverse of the momentum relaxation time. This is 
certainly a very intuitive result; however, let us sketch 
a formal proof for this assertion. The time-dependent 
spin-current operator in the presence of Rashba coupling 
reads 
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where the time evolution includes impurity scattering. 
To lowest order in the spin-orbit coupling and the impu- 
rity scattering we have 



III. ELECTRONS WITH RASHBA COUPLING 

Sinova et al. 15 have considered the spin-Hall effect of 
non-interacting electrons confined to the two-dimensional 
(xy-) plane of a quantum well and being subject to 
Rashba spin-orbit coupling 2 1 ^ 3 , 
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where to* is the effective mass, a the Rashba coefficient, 
and the other notations are standard. We note that in 
systems where both the Rashba and the Dresselhaus spin- 
orbit coupling 22 ' 3 are present, various interesting trans- 
port effects can arise from the interplay of these two 
terms, for recent studies see e.g. 5 ' 7-9 ' 11 ' 12 . For simplicity, 
however, we shall concentrate here on the Rashba term 
only. The Hamiltonian (3) has two energy branches, 
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where the time evolution of <7q is only due to the Hamilto- 
nian (3) and evaluated in the above expression (6), while 
p x (t) contains the impurity scattering but not the spin- 
orbit coupling. Now it is useful to note that, in order 
to compute the expectation values in the Kubo formula 
Eq. (1), only matrix elements of the time-dependent mo- 
mentum operator p x (t) which are diagonal in the wave 
vector index are needed. This enables to apply superop- 
erator techniques developed in Rcfs. 24 yielding 

(p°Mtt « (e- n °*Px(0))^ = (e- t/T p° x (0)) m , (10) 

where fio is the scattering master operator in lowest or- 
der of the scattering potential 24 . It is the same operator 
as it occurs as the scattering term in the usual Boltz- 
mann equation when evaluated in lowest oder via Fermi's 
golden rule. For impurity potentials being isotropic in 
real space, the momentum p x is an exact eigenfunction 
of SIq , and the eigenvalue is given by the well-known 
inverse momentum relaxation time l/r(e) 24 ' 25 which in 
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general depends on the energy s(k). To lowest order in 
the Rashba coupling, this energy argument can be re- 
placed with the Fermi energy in the absence of spin-orbit 
interaction. We note that this momentum relaxation rate 
1/r is the same as obtained in the standard diagram- 
matic approach and thus contains the vertex correction 23 . 
However, this vertex correction vanishes for short-range 
isotropic scatterers. The above argumentation refers to 
the Rashba Hamiltonian (3) for conduction band elec- 
trons; similar considerations can be performed in the case 
of valence band holes to be discussed further below. 

For u> = 0, but finite momentum relaxation rate 1/r > 
0, Eq. (6) yields 



<,*(0) = -<^(0) = - - 



■ tan 



, £ R £ f 
\h/r) 




(11) 



where we have introduced the "Rashba energy" e R = 
m*a 2 /h 2 . Clearly, this is the energy which has to be 
compared with the energy scale %/t of the impurity scat- 
tering in order to obtain the correct zero-frequency limit 
of the spin Hall conductivity. If the impurity scattering 
is weak compared to spin-orbit coupling, ert/Ti 1, we 
have the expansion 
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where we have additionally assumed that the Fermi en- 
ergy £f is larger or at least of the same order as e r , 
which is usually the case in experimental situations. The 
zeroth-order contribution e/Sir is the result obtained by 
Sinova et al. using directly a zero-frequency perturbative 
expression for the spin-Hall current neglecting effects of 
a finite electron quasi-particle lifetime, (cf. Eq. (8) in 
Ref. 15 ). Remarkably, this value is universal in the sense 
that it is independent of a. Therefore, it predicts a finite 
spin-Hall conductivity even in the limit of vanishing spin- 
orbit coupling, a — > 0, which is certainly an unphysical 
feature. However, this paradox can be resolved by the 
observation that the above two limiting processes do not 
commute. In fact, in the opposite limit £rt/Ti <C 1 the 
lowest order of the second term on the r.h.s. of Eq. (6) 
cancels the first one, and the spin-Hall conductivity is 
given in leading order by 
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Thus, to obtain the correct dc spin-Hall conductivity, the 
"Rashba energy" er — m*a 2 /h 2 should be compared 
with the energy scale %/t of the impurity scattering. If 
Eji> H/t the spin-Hall conductivity is close to its "uni- 
versal" value e/87r, while it vanishes for small spin-orbit 
coupling and finite impurity scattering. In epitaxially 



grown GaAs quantum wells mobilities fi = er/m* of or- 
der 100m 2 /Vs can routinely be achieved, corresponding 
to values for H/t of order O.OlmeV. This is safely smaller 
than typical values for the Rashba energy reported from 
experiments 26_34 being of order 0.1 .. . l.OmeV. However, 
it should be noted that the Rashba coefficient is typi- 
cally proportional to an external electric field applied in 
the growth direction of the quantum well. Therefore also 
smaller values of the Rashba energy are possible where 
the finite momentum relaxation time will influence the 
value of the spin-Hall conductivity. 



IV. HOLES IN THE VALENCE BAND OF III-V 
SEMICONDUCTORS 

Murakami, Nagaosa, and Zhang have investigated 
spin-Hall transport in three-dimensional bulk systems 
of holes in the valence band of III-V semiconductors 14 . 
These authors used a phcnomcnological semiclassical the- 
ory to describe adiabatic hole dynamics. Their work 
was revisited most recently by Culcer et al. 16 within the 
framework of a semiclassical theory of wave packet dy- 
namics. Here we will evaluate the spin-Hall conductiv- 
ity using the rigorous Kubo formula (1). Our starting 
point 1 ' 16 is Luttinger's four-band Hamiltonian for heavy 
and light holes in the spherical approximation 35 , 
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Here m is the bare electron mass, and S are spin-3/2- 
operators. The dimensionless Luttinger parameter 71 
and 72 describe the valence band of the specific mate- 
rial with effects of spin-orbit coupling being included in 
72. The eigenstates of (14) can be chosen to be eigen- 
states of the helicity operator A = (k ■ S)/k. The heavy 
holes correspond to A — ±3/2, while the light holes have 
A = ±1/2. From the Kubo formula (1) one finds for the 
frequency-dependent spin-Hall conductivity after lengthy 
but elementary calculations 
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are the Fermi wave numbers for heavy and light holes, 
respectively. Again it is instructive to consider the case 
for weak spin-orbit coupling, 72 <C 71. For 72 = we 

have fcj = k l f =: kp = y ^2ms / / '71ft 2 , and therefore the 
integral in Eq. (15) vanishes for finite 1/r > and all 
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frequencies lu. Thus, as before for the case of electrons, 
the dc spin-Hall conductivity vanishes for vanishing spin- 
orbit coupling if a finite momentum relaxation rate is 
taken into account. This result is in contrast to state- 
ments in Refs. 14 ' 16 , where such dissipation effects were 
neglected. Specifically, for u = we have 
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The remaining integral is elementary leading to a rather 
tedious expression which shall not be given here. How- 
ever, we see that the value of the above integral is gov- 
erned by the ratio of Ti/t and the "spin-orbit energy" 
e so := 7i 2 72(fc°) 2 /m = 25/72/71, since k'j is a typical 
wave number in the integration interval. If Ti/t 3> e so 
the spin-Hall conductivity vanishes as 
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where we have also assumed that the ratio 72/71 is small 
as it is usually the case 36 . In the opposite case Ti/t <C e so 
one finds 
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We note that the zeroth order of this result agrees with 
the expression given in Rcf. 16 for the dc spin-Hall con- 
ductivity neglecting dissipation effects 37 , but differs from 
the result reported in 14 . On the present stage we can- 
not comment on the question whether this difference is 
an artifact of the semiclassical approache used in Ref. 14 , 
whereas the present result is obtained from a rigorous 
linear-response theory given by the Kubo formula with 
full frequency dependence. 

Let us illustrate our results on the typical example 
given in Ref. 14 , where a GaAs sample with hole den- 
sity n = 10 19 cm -3 , corresponding to a Fermi energy of 
order a few ten meV, and a mobility of /i = er/m* = 
50cm 2 /Vs. To obtain an upper bound for r we take m* 
to be the heavy-hole mass, m* w 0.5m, corresponding to 
71 ~ 7 and 72 w 2.5 for GaAs 36 . This leads to a lower es- 
timate for Ti/t being also of order a few ten meV. Thus, 



in the above scenario, the finite momentum relaxation 
time r must be taken into account when calculating the 
spin-Hall conductivity, differently from the approach in 
Rcf. 14 . 



V. CONCLUSIONS 

We have studied the spin-Hall transport of electrons 
and holes in semiconductors using the Kubo formula in 
the correct zero-frequency limit taking into account the 
finite momentum relaxation time of carriers in real semi- 
conductors. This approach allows to analyze the range of 
validity of recent theoretical findings 14 ^ 16 . In particular, 
the spin-Hall conductivity is found to vanishe for vanish- 
ing spin-orbit coupling if the correct zero-frequency limit 
is performed. In the case of conduction band electrons 
in the presence of Rashba spin-orbit coupling in a high- 
mobility quantum well, spin-orbit interaction dominates, 
for typical experimental for the rashba coefficient, the 
effects of momentum relaxation, and the spin-Hall con- 
ductivity is close to its "universal" value as predicted in 
Rcf. 15 . This situation can be different for typical p-doped 
bulk samples, where dissipation can substantially affect 
the spin-Hall transport. 
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